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In the standard simple treatment of magneto-ionic theory the effect of collisions is 
allowed for by assuming that an electron experiences a retarding force proportional to its 
velocity, and the two complex refractive indices of an ionized medium are then given by 
the Appleton-Hartree formula. Experiments have shown, however, that the collision 
frequency is approximately proportional to the square of the electron's velocity, and proper 
allowance for this, using the Boltzmann equation, leads to a modification of the Appleton- 
Hartree formula which was given by Sen and Wyller. In this tutorial paper the same 
modified formula is derived, but by a different method which, it is hoped, can be followed 
by readers not intimately familiar with the previous literature. Some numerical results 
are presented in which the modified and unmodified formulas are compared. It is concluded 
that the standard Appleton-Hartree formula can be used without modification for nearly 
all radio propagation problems in the ionosphere, provided that the correct effective value 
of the collision frequency is used. The modifications may be important, however, in the 
theory of wave interaction and for waves of very low frequency whose wave normals are 
perpendicular to the earth's magnetic field. 

1. Introduction 

The two refractive indices for radio waves in an ionized medium are given by the Appleton- 
Hartree formula [Appleton, 1932; Ratcliffe, 1959]. They are, in general, complex numbers 
whose values are affected by the collisions of the electrons in the medium with other particles. 
In the simplest treatments, this is allowed for by the method of Lorentz [1952, sec. 120 and 
note 57; see also Ratcliffe, 1959, ch. 4, 5] in which it is assumed that an electron of mass m 
experiences a retarding force —mvx proportional to its velocity v, where v is the effective 
collision frequency and is assumed to be a constant. But experiments (Huxley, 1959; Phelps 
and Pack, 1959] show that the average number of collisions v which an electron makes per unit 
time with other particles, depends on the electron velocity v, and when this is allowed for, a 
more elaborate treatment of magneto-ionic theory is needed. The effect of collisions of elec- 
trons with neutral molecules and with charged ions is somewhat different [Ginzburg, 1964]. 
In the lowest part of the ionosphere the concentration of molecules is so great that collisions of 
electrons with molecules dominate, and collisions with charged ions can be ignored, and only 
this case is considered in the present paper. There is good evidence [Huxley, 1959; Phelps 
and Pack, 1959] that the average number of collisions of an electron with molecules per unit 
time is proportional to the square of the velocity. 

This problem has already been discussed by Molmud [1959] and more fully discussed by 
Sen and Wyller [1960] whose results are the same as those given here. Thus there is nothing 
new in the formulas derived in the present paper, but the method of derivation is different. 
Sen and Wyller used a coordinate system with oblique axes, and they allowed for the harmonic 
time variation of the radio wave fields by including real factors cos ut and sin u>t. In the present 
paper we use orthogonal axes and a complex time factor exp (ico t) as was done, for example, 
by Allis [1956], and these lead to considerable simplification. 

The results are important in studying the propagation of radio waves in the lowest part of 
the ionosphere. In one kind of experiment Gardner and Pawsey [1953] have observed the 
reflection of vertically incident radio waves from small discontinuities of electron density in 
the height range 65 km to 90 km. By measuring the ratio of the reflected amplitudes for the 
ordinary and extraordinary rays it was possible to deduce the effective collision frequency. 
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But Belrose and Burke [1964] have shown that the results are significantly different for the 
simplified Lorentz theory and for the more exact theory, which leads to the formulas of Sen 
and Wyller. The more exact theory may also be important when studying the reflection of 
waves of very low frequency from the ionosphere, for in some published work [Budden, 1955; 
Barron and Budden 1959] the simple Lorentz treatment was used. The modifications intro- 
duced by the more exact theory have been discussed by Johler and Harper [1962] and by Wait 
[1962, appendix B]. But recently Deeks (private communication) has calculated the reflection 
coefficients of the ionosphere for a range of very low frequencies, using both the modified and 
unmodified formulas, and finds that the differences are negligible (see sec. 9 below). 

The treatment given both by Sen and Wyller and in the present paper is based on the 
Boltzmann equation, which is set out in section 3. It is simplified by assuming that the field 
of the radio wave is very small and gives only a small perturbation j x of the distribution func- 
tion. This leads to a linear equation for ^ given at the end of section 3. It contains a term 
which allows for collisions, and is discussed in section 4 where it is expressed as an integral — 
the collision integral. The equation for/i is then solved in section 5, and the solution is used 
in section 6 to derive the electric susceptibility matrix of the medium. From this the refrac- 
tive indices can be derived, and this is done first for an isotropic medium in section 7. The 
application to the full magneto -ionic theory is discussed in section 8 where Sen and Wyller's 
modifications to the refractive indices and wave polarizations are derived. Some numerical 
results are presented and discussed in section 9. When the field of the radio wave is stronger, 
some of the terms which are neglected in the simple theory are no longer small, and the con- 
sequences of this are very briefly reviewed in section 10. 

2. Notation 

The formulas in this paper are valid in any system of self-consistent rationalized units. 
The principal symbols are as follows: 

B= magnetic induction of superimposed constant magnetic field. 

c = velocity of electromagnetic waves in free space. 
E= electric field of radio wave, components E x , E y , E z . 

e = charge of electron. 

j= Boltzmann distribution function for electrons. 
H= magnetic field of radio wave, components H x , H y , H z . 

J = current density. 
ifr=Boltzmann's constant. 

Jc = u/c, 2t divided by wavelength in free space. 
m=mass of electron. 

N= number of electrons per unit volume. 
N m = number of molecules per unit volume. 
n= refractive index. 

n= vector of length n in direction of wave normal. 
P= J/ico, electric polarization. 

reposition vector, components x, y, z. 
d 3 r=dxdydz. 
T= absolute temperature of electrons. 

/=time. 
U=l-iZ. 

L/ e ff=l ^eff- 

v= electron velocity vector, length v, components v x , v y , v z . 
d*\=dv x dvydv z . 
X=Ne 2 /e Q mu 2 . 

Y=eB/com, vector of length Y. 
Z=v(v)/ca. 
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Z m =v m /(a=l/w. 

5 
Z eff =-v TO /w, effective value of Z used in standard magneto-ionic theory. 

Zq= characteristic impedance of free space. 
e=scalar dielectric constant. 

c=tensor dielectric constant, with elements e u , (i, j=x, y, z). 
e h e2j e 3 =principal axis values of dielectric constant; e 3 refers to axis parallel to B. 
6= angle between vector Y and wave normal. 
v = collision frequency, a function of electron velocity. 

2 
p m =z-X average value of v, when i>ocz? 2 . 

o 

Pi p 0? p x =wave polarizations in magneto-ionic theory. 
co = 2 7rX frequency of wave. 
grad r = vector operator in ordinary space, with components d/dx, d/d//, d/ds. 
grad v = vector operator in velocity space, with components d/c)o x , b/5?' /y , d/dtf*. 

3. Linearized Boltzmann Equation 

We consider a neutral, ionized, homogeneous medium containing N electrons and N 
positive ions per unit volume, and with neutral molecules whose concentration N m is very 
large compared with A" so that electron-electron and electron-ion encounters can be ignored. 
A molecule is very massive compared with an electron so that as a first approximation the 
molecules can be assumed to be infinitely massive and at rest. Thus the state of the medium 
can be described by the Boltzmann distribution function /(r, v) for the electrons. There is 
a superimposed constant magnetic field of induction B. The function / must satisfy the 
Boltzmann equation 

g+vgrad r /+agrad T /=^ (1) 

[see, for example, Chapman and Cowling 1958, ch. 3] where a is the acceleration imparted 
to an electron by all forces except encounters. The term -~- is assumed to arise only from 
electron-molecule encounters. 



The acceleration is given by 



a=^-{E+vAB} (2) 



where E is the electric field of the radio wave, which is assumed to be so weak that the force 
exerted on an electron by the magnetic field H of the wave can be neglected. It is now further 
assumed that when E is zero, and when the electrons have reached a state of equilibrium, 
the function/ is given by the Maxwellian distribution function 

fo(v)=N (^EfT exp (- w " 2 /2#r). (3) 

This is independent of time so that both -^r and -4~ are zero. In the ionosphere the electron 

ot of 

concentration N in (3) is a function of the height z and might give a contribution to the term 
v-grad,. /' in (1). But it can be shown that the order of magnitude of this term is small enough 
to be neglected for frequencies greater than about 1 kc/s, though it could possibly become 
important in the ELF range (less than about 1 kc/s). Thus we shall assume that N is in- 
dependent of r so that the medium is homogeneous. 

It follows from (3) that grad v /ois a vector parallel to v. Thus its scalar product with 
the term vAB in (2) is zero. The field E is assumed to be very small, and when it is present 
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the function/ undergoes a small change so that we now have 

/(T,F,*)=/o+/l(V,T,0 (4) 

where the perturbation/! is so small that squares and products of/ and E can be neglected. 
If these results are used in (1) we obtain 

^+v ■ grad r .A+^ E grad v /„+^ (yAB) grad,/,^- (5) 

This equation is linear in E and/. The electric field E and magnetic field H of the radio 
wave must satisfy Maxwell's equations which are also linear. We seek a solution of the 
equations which represents a monochromatic wave of angular frequency co. It follows from 
the linearity of the equations that a solution exists in which/, E, and H all vary harmonically 
in time with angular frequency w. This may be allowed for by assuming that each of them 
includes a complex time factor exp (i&t) which is customarily omitted from the equations, 
and when interpreting the complex f h E and H it is understood that the real parts are to be 
taken. Thus in (5) we may set 

g=W, (6) 

The term -~i is considered in the next section. 
ot 

In magneto-ionic theory we seek a solution which represents a plane monochromatic 

wave. Thus /, E, and H all contain a factor 

exp i(ut— &n • r) (7) 

where the vector n has the direction of the wave normal and (complex) magnitude equal to 
the refractive index. If this is used in the second term of (5) it gives 

v • grad r /i=— ikn • yf lm (8) 

This term depends on how/i varies in space, and its effect on the results is sometimes called 
"spatial dispersion" [Ginzburg, 1964]. It is very small compared with (6) provided that 

n«fi/vi (9) 

Now in magneto-ionic theory we are usually concerned with refractive indices of the order of 
unity or less, whereas the number of electrons whose velocity v is comparable with c is very 
small. Thus (9) is satisfied and we may neglect the spatial dispersion term (8). By retaining 
the term (8) it is possible to study plasma waves in the medium, but this is not part of classical 
magneto-ionic theory. The term (8) also leads to Landau damping which is negligibly small 
for radio waves when n is of order unity. This term has a very small effect on the polarizations 
and refractive indices of magneto-ionic theory [see Ginzburg 1964, sec. 12]. In the rest of 
this paper it is neglected, except for a brief mention at the end of section 5. 
The third term in (5) can be evaluated from (3) and gives 

^ E grad v fo=-j£j, E . v/ . (10) 

With these simplifications (5) gives 

wtfi-gy E • v/„+^ (TAB) • grad v /i=^ (11) 

This is the form of the linearized Boltzmann equation which will be used in the rest of this paper. 
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4. Collision Integral 

The Boltzmann distribution function /(r, v) is defined so that 

/(r, v)r/ : "r(/ 3 v 



(12) 



is the number of electrons in the small volume d h r of ordinary space, having velocities within 
the small cell d s \ of velocity space. Thus the electrons of this small group all have a velocity 
very close to a given value v. When one electron of the group makes a collision it is removed 
from the cell d 3 \. Other electrons in different cells make collisions which sometimes result in 
an electron moving into the cell d z v. It is these two processes which contribute to the term 

-~rr and an expression for this must now be derived. 

It is assumed that a collision is practically instantaneous so that E and /^ do not change 
appreciably while it is occurring. Thus the duration of a collision must be very small compared 
with 1/co. This is believed to be true for electron-molecule encounters when co refers to a radio 
wave of low or medium frequency. It may not be true for electron-ion encounters or for very 
high frequencies and microwaves. The problem is then more difficult and it is not always 
possible to use the Boltzmann equation in these cases. 

The molecules are assumed to be infinitely massive and at rest. Thus an electron which 
approaches a molecule is deflected by the encounter but the magnitude of its velocity is un- 
changed. It jumps from the cell d 3 v to another cell d 2 \ f on the same sphere, of radius v, in veloc- 
ity space. Suppose that there is a molecule at the point M in ordinary space (fig. 1), and t hat 
an electron is approaching it along the line AB parallel to the velocity v. The volume of the 
cell cZ 3 r in ordinary space, though small, is very large compared with the region represented by 
figure 1 . This cell thus contains a large number N m d 3 r of molecules, where N m is the number 
of molecules per unit volume. An electron which makes an encounter with one of these mole- 
cules remains in the cell dh in ordinary space, though it jumps to a different cell in velocity 
space. Now imagine a plane drawn through M normal to AB, and take polar coordinates 
b, a in this plane with origin at M, radius b and angle a measured from some reference line MC. 
The radius b is called the impact parameter. Then PQRS is a small area bdbda of this plane. 
The probability that an encounter occurs in time dt, with an electron approaching the area 
PQRS, is equal to the probability that there is an electron in the volume bdbdavdt of ordinary 
space, so that the probability is 

J{y)d 3 \bdbdavdt. 



Figure 1. Direct encounter of an electron with velocity 
v approaching an infinitely massive molecule at rest 
at M. 

The thick line represents the path of the electron. 
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After the encounter this electron has a new velocity v' where v'=v but v'=^ v. 
of this kind of encounter for all the molecules in the volume cPr is 



The probability 



N m yf(y)d?\d*rbdbd<xdt. 



(13) 



Electrons also enter the cell d?\ through collisions which are called "inverse encounters." 
There will clearly be some electrons which, after an encounter, are traveling with velocity v 
so that they are then in the cell dH. Some examples are illustrated in figure 2b, in which the 
electrons initially have a velocity \' or v" or v"' all with the same magnitude v but with direc- 
tions which may be anywhere within the solid angle Aw. Of these we consider only those 
electrons whose final velocity vector v, when projected backwards, intersects the area P'Q'R'S' 
of figure 1, which lies near the point with polar coordinates b, a + 7r. Thus for these collisions, 
the impact parameter b is the same as for the direct encounter considered in (13), and the angle 
through which the electron is deflected is the same for both cases. Hence the initial velocity 
vector v' of the inverse encounter must be the same as the final velocity vector v' for the direct 
encounter, and the electron orbits for the two encounters are identical in shape, though differ- 
ently oriented in space. The process of counting these inverse encounters is exactly the same 
as counting the direct encounters, which led to (13). If a small change is made in v so that it 
remains within the cell d 3 v of velocity space, then there is a resulting change in v' which remains 
within a corresponding cell c/V. Now the volumes d 6 \ and c/ 3 v' of the two cells are equal. 
This can be seen by expressing v and v' in spherical polar coordinates, v,&,<p and v §v &',<p' respec- 
tively, using as polar axis the bisector of the angle between v and v'. Then v=v', &=&' , <p= 
<p' + t whence d?\=v 2 sin &dvd&d<p=d?\' . This is a special case of the principle of detailed 
balancing and a formal proof, for a more general type of binary collision, is given by Chapman 
and Cowling [1958, sec. 3.52]. 

The number of electrons leaving the cell d s rd d v in time dt is given by (13). The number 
entering the cell from the inverse encounters is the same expression with /(v)cPv replaced by 
f(v')dW =j(y')d?v. The resultant number entering d 3 vd s r in time dt is thus 



N n {f(y')-f(y) }vd z \d?vbdbdadt. 



(14) 



K1 



To find -^j-t equation (14) is divided by d 3 \d 3 rdt and then integrated with respect to b 
and a. Thus 

M 

dt 



d=N m vf*"£lf(v')-f(Y) }bdbda. 



(15) 



This is the collision integral. In the integrand the vector v is a constant, and v' and v have 
equal magnitudes. But the direction of v' is a function of b and a. 



Figure 2. In (a) a direct encounter is shown, in which 
an electron initially with velocity v leaves the cell d 3 v 
centered on the point v in velocity space. 

Some inverse encounters are shown in (b), in which electrons with 
various initial velocities v', v", v'" all enter the same cell d 3 v after 
the encounter. 




(a) 
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When the impact parameter b is large, it usually happens that the electron is almost 
undeflected by the encounter, so that v'«v and the integrand is very small. This occurs for 
electron-molecule encounters, where the forces are of short range, and the integral (15) con- 
verges as b increases. Thus the upper limit for b is taken as infinity. But for electron-ion 
encounters the inverse square law of force is of longer range and the integral does not converge. 
Then the upper limit for b must be assessed by considering the screening effect of the neighbor- 
ing plasma. In the present paper, however, these considerations need not concern us since, 
following Sen and Wvller [1960], we shall use experimental results to give the final form of the 

term -^7- 
ot 

When /(v) is independent of the direction of v, the integrand of (15) is zero. Thus / (v) 

contributes nothing to the collision integral, and the/ in (15) may be replaced by f lm Equations 

(11) and (15) then give 

ico/i-^E.v/ +-(vAB) .grad v f x =N m v P* P{ /i(v / )-/ 1 (v) }bdbda (16) 

i±± in> Jo Jo 

winch applies for each value of the velocity vector v, and is to be solved to find the function/i(v). 

5. Solution of the Equation 

Equation (16) must apply at each point v in velocity space and we shall study it for those 
points which lie on a sphere of fixed radius v. Then/j is a function of the direction of v, and 
may be expanded in a series of spherical harmonics thus 

/i(v)=/ioW+/n(v)+/i 2 (v)+ . • • (17) 

where /^(v) is a spherical harmonic of order r. This is substituted into (16) and the terms 
representing harmonics of the various orders are equated. The second term in (16) is a spherical 

harmonic of the first order, which could be written —j==,jj£v cos $ where $ is the angle between 

E and v. The polar axis is then the direction of E and may be a complex direction since E in 
general represents an elliptically polarized wave [Budden, 1964b]. 

Now the term f 10 (v) is independent of the direction of v, so that grad v /i (fl) is a vector 
parallel to v and thus contributes nothing to the third term of (16). Similarly /io(fl) contributes 
nothing to the collision integral because fio(v')=jio(v). The only spherical harmonic of zero 
order in (16) is in the term ico/i, so that 

/io(t»=0. (18) 

Next consider the spherical harmonic of first order / n (v) which proves to be the most 
important term. It may be written 

/„(T)=V.V(») (19) 

where the vector V depends only on the magnitude of v and not on its direction. Now 

grad v (v.V) = V+(vg)vA>. (20) 

The last term of (20) is a vector parallel to v and contributes nothing to the third term of (16), 
which therefore gives 

£ (TAB).grad T /„=£ (v A B)-V=^ (BaV)-t (21) 

which is a spherical harmonic of the first order. The contribution of (19) to the collision 
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integral is 



N m vY- (^ ( W (v'-v)bdbda. (22) 

Jo Jo 



To evaluate this, choose Cartesian axes v X) v y , v z in velocity space so that v x is in the reference 
direction MC in figure 1, and v z is in the direction of v. Suppose that in the encounter the 
electron is deflected through an angle x(b). Then the Cartesian components of v' — v are 

— v sin X cos a,—v sin X sin a, v (cos X— 1). (23) 

Integration of the first two with respect to a gives zero so that only the v z component remains 
and this is parallel to v. Hence (22) gives 

-2wN m v\'Y f°°(l-cos x)bdb. (24) 

Jo 

A(v)=2w f °°(1— cos X)bdb. (25) 



Now let 



This is called the "transport area of cross section" for the encounter. It may depend on v, 
but is independent of the direction of v. It is determined by the mechanics of an electron- 
molecule encounter. For example, it is easy to show that if the molecules are rigid spheres of 
radius a, then A = wa 2 . 

If each molecule presents an effective area A(v) to the oncoming electrons, then the 
probability of an electron encountering a molecule within a path dx is N m Adx. Thus the mean 
free path is 

\(v) = l/{N m A(v)} (26) 

and the average number of collisions per unit time is 

v(v)=v/\=N m Av. (27) 

Equation (27) is adopted as the definition of the collision frequency v(v). Huxley [1959] has 
described experiments which show that A is proportional to v, so that v(v) in (27) is proportional 
to v 2 , and this result is used later. If (27) is combined with (25), the contribution (24) to the 
collision integral is 

—\'Yv(v) (28) 

which is a spherical harmonic of the first order in velocity space. It is shown below that 
the spherical harmonics of second and higher order in (17) do not contribute any terms of 
first order in (16). Hence the spherical harmonics of first order give 

io.V-T-^/ E-v+^(BAV)v+,V-T=0 (29) 

which can now be solved to find V. This is done in the following section. 

The spherical harmonic of the second order in (17) may conveniently be written in tensor 
notation thus 

f l2 (Y)=W ij (v)v i v j (i,j=x, y, z) (30) 

where W tj is a symmetric tensor whose trace is zero, and the summation convention for repeated 
suffixes is used. The contribution of (30) to the collision integral in (16) is 

N m vWij [(v'tVj—ViV,) bdbda (31) 

and here the integral can be evaluated by using the same Cartesian coordinate system as was 
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used above to evaluate (22). The elements of the integrand of (31) can then be written in terms 
of the polar coordinates v, x, & of v ' (compare (23)) and integrated with respect to a. In this 
way it can be shown that (31) gives 



-NnM 



r°° 3 

T ijViVj o sm2 x bdb (32) 

Jo ^ 



which is a spherical harmonic of the second order and contributes nothing to the first order 
terms in (29). 

The contribution of (30) to the third term in (16) may similarly be found, since on converting 
to tensor notation 

a k = {grad v / 12 (v)U=^- {W^v^} =2W tk v ( + ViV f k *£« (33) 

and 

(vAB) •a=e m B J a k v i (34) 

where e ijk is the third order isotropic tensor [Jeffreys and Jeffreys, 1956, ch. 3]. Substitution 
of (34) into (33) gives for the third term of (16) 



-j 2€ m B J W lk v i v l +e ijk B j -^v i v k v l vJv V (35) 



in which the first term is a spherical harmonic of order two, and the second term is of order 
four, and there is no term of the first order. 

More generally it can be shown that a spherical harmonic of order /' in (17) does not 
contribute terms of smaller order in the third term of (16). By expressing it as Legendre 
functions and associated Legendre functions it can further be shown [Allis, 1956, p. 409] that 
it contributes only a spherical harmonic of the same order in the collision integral. Thus the 
spherical harmonic terms of order r=2 or greater, in (17), when inserted into (16) give a set of 
linear homogeneous equations which do not contain the field E. Since we are only interested 
in the disturbances produced by the radio wave field E, we may clearly take the terms /i r (v) 
to be zero when r>2. 

It may be noted that this simplification does not occur if the spatial dispersion term (8) 
is included, for then the spherical harmonic of first order, /n(v), when inserted in (8), gives a 
spherical harmonic of order two. Since / n (v) = V-v depends on E, as is shown by (29), this 
introduces a term depending on E into the second order spherical harmonic component of (16). 
Thus it is no longer permissible to take IT ?7 = 0, which leads to some complication, but it is 
beyond the scope of this paper. 

6. Susceptibility Matrix of the Medium 

In the last section it was shown that only the second term of (17) is to be retained and 
it is given by (19). Thus the distribution function (4) for electrons becomes 

./=/o+V-v (36) 

where the vector V(v) is given by (29). Every term of (29) is the scalar product of v with some 
other vector, and since the equation must be true for all v, we may omit the factor • v and obtain 

{^+K*>)}V+^BAV=-^ f ./oE. (37) 

Now it is convenient to use the standard notation (U. R.S.I. 1956) of magneto-ionic theory, 
in which the vector B is denoted by Y and is antiparallel to B because e is a negative number. 
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Similarly v/oi is denoted by Z. This is slightly different from the usual Z because here v{v) 
and Z{v) may depend on the electron velocity. Finally we use U(v) = l— iZ. Then (37) gives 



l7V— rt7\V= 



iuKT 



/oE. 



(38) 



This is very similar to an equation used in standard magneto-ionic theory [Budden, 1961, 
eq. (3.15)] and may be solved in the same way. Choose Cartesian axes with the 2-axis parallel 
to Y and write (38) in matrix form thus 



or 



T x 



-iY u o 

v.0 Us 

By inversion of the 3X3 matrix this gives 



<V*s 



e 



r E x "\ 

E 



\E*s 






i 



- U ( U 2_ Y 2) 



iUY U 2 



(i 



» 





U 2 -Y 2 






iuKT 



E y 

Ezj 



(39) 



(40) 



This is more concisely written by using the matrix notation 

<*U 2 -iUY 



Q 



l 



U{U 2 -Y 2 ) 



0^ 





so that (40) gives 



iUY U 2 
L0 U 2 -Y 2 J 



V =A QE - 



(41) 



(42) 



The current density vector J in the ionized medium is given in terms of the Boltzmann 
distribution function thus 



[=e JIf /(V)V ^ 



(43) 



the integral being taken through the whole of velocity space. To this integral the symmetric 
function / (v) contributes nothing so that insertion of (36) gives 



J = e 



JJI <vV) 



v-dH. 



(44) 



Now use spherical polar coordinates v, &, <p in velocity space and choose the direction of V as 
the polar axis. Then the components of (44) perpendicular to V are easily shown to be zero 
so that (44) gives 



J= e 



(*2tt fir f°° A r°° 

\(v)i A cos 2 # sin #dvdU<p=- ire \(v) 
Jo Jo Jo w Jo 



•vSlv. 



(45) 



Instead of J it is convenient to define an electric polarization vector 
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(46) 



as is done in standard magneto-ionic theory [Ratcliffe, 1959; Budden, 1961]. Then a combina- 
tion of (46), (45), and (42) gives 

P=e M.E, (47) 

where 



M 



is called the susceptibility matrix of the medium. 

In standard magneto-ionic theory it is assumed that v and Z are independent of the electron 
velocity v, and then Q is independent of v and may be taken outside the integral in (48). Now 
Jo is given by (3) so that the integral may be evaluated and the result is 

M=-ZQ (49) 

which is the same as in standard magneto-ionic theory [Budden, 1961, eq (3.24)] and leads to 
the Apple ton-Hartree formula. 

7. Isotropic Ionosphere 

For an isotropic ionosphere the earth's magnetic induction B is neglected so that Y is 
zero and Q in (41) becomes a scalar 1/U. Then the susceptibility M also becomes a scalar M 
which is related to the complex dielectric constant e, and the complex refractive index n thus: 



n 2 =e=l+M=l-- 



4tt6 2 r°Mvw 



r^*- (5o) 



3€ O>2KT, 

Again if U is independent of v this leads to the standard result 

n 2 =l-X/U. (51) 

Now it might be expected that when an electron moves in a gas of infinitely heavy molecules 
which are at rest, the probability of its encountering a molecule would depend simply on the 
position and direction of its path, and not on its speed. Thus the mean free path would be 
independent of velocity, but the average number of collisions made in unit time would be pro- 
portional to the velocity. This case has been considered by Margenau [1946] and by Ginzburg 
[1964, sec. 6] who gives curves showing how the real and imaginary parts of e depend upon 
frequency co. 

The true position is less simple because fast moving electrons can more easily excite the 
rotational states of molecules, so that significant collisions occur more easily for fast electrons 
[Huxley, 1959]. Careful experiments [Crompton, Huxley, and Sutton, 1953; Huxley, 1959; 
Phelps and Pack, 1959] have shown that the collision frequency v(y) is more closely propor- 
tional to the square of the velocity. Hence following Sen and Wyller [1960] it will be assumed 
that 

V {v) = ^vK (52) 

3 
By using (3) it can be shown that the average value of v is then - v m . Now let 

u 2 =v/v m =mv 2 /2KT. (53) 

By using this with (3) and (50): 

, , AX C" u 3/2 e-"du .... 

n 2 =«=l— 5-7= = : ;-■ (54) 

3Vt Jo l—lUv m /u 
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This integral may be expressed in terms of standard integrals defined and tabulated by 
Dingle, Arndt, and Koy [1956-1957], namely 

#! Jo ^ +w 2 

Now (54) is separated into its real and imaginary parts, and with 

w=u/v m (56) 

we obtain 

n 2 =€ ^ l+M=1 _ X {w 2 Gz / 2(w)+liwe5 / 2(w)} (57) 

which replaces the standard formula (51). 

When vju is very small the integral in (54) may be evaluated approximately by using 
a binomial expansion for (1— iuv m loi)~ l and integrating term by term. This gives 



■»-*{i-t(?)+-}-*65+-> 



(58) 



When the formula (51) of standard theory is used, it is assumed that Z can be given an 
effective value Z eff which is independent of electron velocity. Then (51) gives 

n 2 =6 = l-X/(l+Z e 2 ff )-iXZ eff /(l + Z e 2 ff ). (59) 

If both v m /oj and Z eff are so small that squares and higher powers can be neglected, (59) agrees 

with (58) provided that 

5 5 

v ett /a)=Z e t =o vj(a=5/(2w), v e u=-^v m ' (60) 

Thus, for very small values of Z the standard formula (59) is approximately correct provided 
that Z and v have the effective values given by (60). Note that v eft is not the average value 
\v m of v, but is f of this average value. For greater values of Z, (59) is in error and (57) 
must be used. To obtain (57) from (59) the real and imaginary parts of e— 1 must be multi- 
plied by correcting factors K r and K u respectively, so that, with (60) 

e=l-7T r X/{l+25/(4w 2 ) }-5iK i /{2w+25l(2w) } (61) 

and by equating this to (57) 

Kr=(w 2 +^ ) (?.«(«); i^=(w 2 +f ) e^(wy (62) 

These two functions are shown in figure 3 plotted with logi (1/Z eff )=logi (fw) as abscissa. 
Similar curves were given by Ginzburg [1964, sec. 6] for the case where v(v) is proportional 
to v. 

When j> m /co=fZ eff is very large (w very small) the integral in (54) may again be evaluated 
approximately by using a binomial expansion for the factor {l+iw/iuvm)}* 1 and integrating 
term by term. This gives 

n*=€~l-! Xw 2 -i \ Xw (63) 

o o 

and by comparison with (61), for small w 

Zr-y; Z,~§. (64) 

The curves of figure 3 approach these limiting values when Z ett is large. 
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Figure 3. The correction factors K r and Kj used for 
converting the real and imaginary parts respectively 
of e — 1, as computed with the simple formula (59), 
to the more accurate values given by (57) . 




8. Refractive Indices of Anisotropic Ionosphere: Modified Appleton-Hartree 

Formula 

For an anisotropic ionosphere Y is not zero and the elements of the matrix (48) must be 
found. Equation (41) shows that three different integrals are involved in (48) with factors 
U 2 , UY, and U 2 — Y 2 in the integrand. Further, the element M zz is identical with Min (50j, 
and is given by (57). The remaining nonzero elements are M xx =M yy and M xy =—M yx , and 
(48) and (41) together show that 



M xx TiM x% 



4?T6 2 r 

~3e a> 2 KTj 



U±Y 



dv. 



(65) 



The two expressions in (65) together with M Z2 comprise the three "principal axis' y values of the 
susceptibility matrix (see Budden, 1961, sec. 3.11, and the references given there; 1964a, 
sec. 2.8), and the corresponding "principal axis" values of the dielectric constant matrix are 



e x =l-\-M xx -iM xyj e 2 = l+M xx +iM xy , e 3 = l+M 22 . 



(66) 



Now assume that v(v) is given by (52) and use the variable u defined at (53). Then, with (3), 
equations (66) and (65) give 



i\_ 4Z p u z '' 2 e~ u i 

2 j 3-v^Jo i±f-mi 



*2 



du 



(67) 



and e 3 is given by (54). The real and imaginary parts of the integrals in (67) may be expressed 
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in terms of w, in (56), and of the Q integrals given by (55). Thus it is easily shown that 

€2 =l-X^(l-Y)w^esf2{w(l-Y)}+^iwe 5/2 {w(l-Y)}'j 
3=1 -X [w 2 CwW+l wgwW] (68) 



€ 3 = 



[note that these e's are not the same as those used by Sen and Wyller, I960]. These three 
complex numbers, involving six G integrals, are sufficient to determine the polarization and 
refractive index for a plane wave. The corresponding values according to standard magneto- 
ionic theory are 

e^l-XKUeu+Y)^ 



e 2 =l-X/(U ett ~Y) 
e,= l-X/U ett 



(69) 



[Budden, 1964a, sec. 2.8] where U et t=l—iZ ett is a constant independent of v. 

The numbers in (68) determine the constitutive relation for the medium, which may be 
written 

e E+P=D=€oeE=€o(l + M)E. (70) 

where e is the matrix dielectric constant. With the coordinate axes used in section 6, in which 
the 2-axis was parallel to the vector Y, the matrix e is given from (66), by 



1(€i+«*) ii(ei-e 2 ) 

—hi{e l —e 2 ) J(«i+«a) 
.0 eu 



(71) 



We now seek a solution of Maxwell's equations which represents a plane wave whose wave 
normal makes an angle 6 with the vector Y. There are several ways of doing this, used in 
standard works on magneto-ionic theory [Ratcliffe, 1959 ; Budden, 1961, 1964a; Sen and Wyller, 
I960]. The following derivation is given only in outline, and uses the author's previously 
published method [Budden, 1961, 1964a]. 

Choose a Cartesian coordinate system in which the 2-axis is the wave normal, and the x-axis 
is in the plane containing the wave normal and the vector Y. In this system the dielectric 
constant e is given by 



f i(*i + e 2 ) cos 2 + e 3 sin 2 0, ^i(e 1 — e 2 ) cos 0, §(€i+€ 2 ) — * 3 sin d cos 6 ^ 



-^i(e l — e 2 ) cos 0, 



«€! + €,), 



K(«i— €2) sin 6 



(72) 



^i(ei + e 2 )— e 3 sin 6 cos 0, i^(ei — € 2 ) sin 0. € 3 cos 2 + i(ei + e 2 ) sin 2 
For the plane wave, the differential operators with respect to the coordinates are 







-=—ikn, 



(73) 



bx by dz 

and when these are used in Maxwell's third and fourth equations, they give, respectively 

nE t =-Z a H x , nE x =Z a H y , H z =0 (74) 
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and 

7lZ Q H y = e xx E x + e xy E y + e xz E z 

-nZ Q H x =e yx E x +e yy E y +e yz E z }- (75) 

e zx E x +e zy E v J i-€ zz E z =0 

The last equation (75) gives an expression for E z which is substituted in the remaining equations. 
The ratio of the two equations in (74) and of the first two equations in (75) then gives 

e yz^zx \ zp \( € yz*zy \ j^ 

*ZZ / \ *ZZ / (n n\ 

P ~W W~7 \ 7 \ — v®' 

■^x tt y I txz^zx \ rp I / *xz*zy \ tp 

so that 

tzz^yx tyz^zxi x^yy^zz ^yz^zy) P 

^xx^zz £xz e zx I \€xy e zz ^xz^zy) P 



(77) 



The ratio p is a complex number which determines the polarization of the wave. Equation 
(77) is a quadratic equation for p and shows that the wave must have one of two possible 
polarizations. If the elements of e from (72) are inserted into (77), then after rearrangement: 

cos0e 3 (e 1 — e 2 ) 

This equation must hold both for standard magneto-ionic theory which uses (69), and for a 
velocity dependent collision frequency when (68) is used. If (69) is used, (78) may be simplified 
and leads to the well-known polarization equation of magneto-ionic theory. In both cases 
the product of the two roots is unity so that 

PoPx=l (79) 

where the subscripts o, x, denote "ordinary" and "extraordinary". This well-known result 
of magneto-ionic theory is therefore unaffected by a velocity dependent collision frequency. 

Next H v from the second equation (74) is substituted in the first equation (75), and E z is 
eliminated by using the third equation (75). If the resulting equation is divided by E x , we 
obtain 



2 K^zz^xx ^xz^zx) I Pytzz^xy ^xz^zy) 



(so) 



and on using the elements of c from (72), this gives 

2= e s( e i + ^)+ip CQS M e i— €2) , 

(e!+€ 2 ) sin 2 0+2 €3 cos 2 6 (<Slj 

Thus the refractive indices n are found by substituting the solutions of (78) for p into (81), 
and this leads to the formula used by Sen and Wyller, namely 

^ e^sin^+^s^ + e^q+cos 2 ^^ 
n (e 1 + e 2 )sin 2 (9+263COS 2 6> & Z) 

This gives the Appleton-Hartree formula if (69) is now used, but it gives the more accurate 
results of Sen and Wyller if (68) is used. Further, n 2 satisfies 

n i {h{e l + e 2 ) sin 2 6+e 3 cos 2 6 }-n 2 { c lf2 sin 2 6+hz(e 1 +e 2 )(l+cos 2 0)} + €l e 2 e 3 =O (83) 

205 



which shows that n is zero when any one of ei, e 2 , or e 3 is zero. Thus in both the Appleton- 
Hartree and the Sen and Wyller case the zeros of n are independent of the direction 6 of the 
wave normal. If v m is not zero, it is not possible for e u e 2 , or e 3 to be zero for any real value of 
X, but the zeros occur for complex values of X and are important when studying the contour 
integrals in the complex A-plane, as used in the phase integral method [Budden, 1961, ch. 20; 
1964a, ch. 5]. 

9. Results 

Many curves were given by Sen and Wyller showing how the real and imaginary parts 
of the refractive indices n depend upon X, which is proportional to the electron concentration N. 
They applied only to the cases 0=0 and 9= %w, that is to purely longitudinal and purely trans- 
verse propagation, respectively, and to the single value Y=2.28. For the ionosphere, in 
which the electron gyrofrequency is about 1.0 to 1.2 Mc/s, this means that their results applied 
only to a frequency of about 500 kc/s. In two of their curves they compared the results of 
their accurate formula using Z m =0A, with those of the Appleton-Hartree formula with Z Gtt = 
|Z m =0.25, as suggested by (60). This would correspond to a value *> W ~3X10 5 sec -1 which is 
about right for a height of 90 km in the ionosphere. In this case the differences between the 
two formulas were significant but small. In other curves they used larger values of Z m namely 
0.5 and 2, and here the differences between the formulas were more marked, although in these 
cases they used Z eU =Z m or \Z m but not |Z m . 

The curves shown here were computed on E.D.S.A.C.2, the digital computer in the 
University Mathematical Laboratory, Cambridge. The two complex values of n and p were 
found for each set of values of A, F, Z. m , and 6 according to the Sen and Wyller formulas, and 
according to the formulas of standard magneto-ionic theory using Z eff =5Z m /2. The C integrals 
in (68) were computed from approximate formulas, given by Hara [1963], which express the 
integrals as rational functions with an error less than 0.7 percent. 

The program was first checked by recomputing the results given by Sen and Wyller 
[I960]. For longitudinal propagation, 6 = 0, the agreement was found to be as good as could 
be judged from the printed curves in their paper. For transverse propagation, 9=ir/2, the 
agreement was good for the extraordinary wave, but for the ordinary wave there was serious 
disagreement. For this case the refractive index should be the same as for an isotropic 
medium as given in section 7. The E.D.S.A.C. program was checked by hand computing a 
few values, and was found to be correct. 

A series of calculations was then made using 0=30°, so that the results would apply for 
radio waves traveling vertically upwards or downwards in the ionosphere in temperate 
latitudes. The first set of values used was F=l/2, Z m =0.02, Z en =2.5 Z m =0.05. They thus 
apply to a frequency of about 2.5 Mc/s and a value j> eff =7.5X10 5 sec -1 , that is to a height of 
about 90 km. The results showed that the differences between the two formulas were too 
small to show clearly on curves of the kind used by Sen and Wyller. A few typical values 
are given in table 1. These suggest that for waves of frequency greater than about 2 Mc/s 
and for heights greater than 90 km in the ionosphere, the standard Appleton-Hartree formula 
may safely be used with v eft =5u m /2. 

Table 1. Values of polarization p and refractive index n as calculated by the Sen and Wyller formula (SW) and 
by the Appleton-Hartree formula (AH), for Y=l/2, 9 = 80°, Z m = 0.02, Z eif = 0.05 







P (SW) 


P (AH) 


rc(SW) 


n (AH) 


X 


Wave 


R 


I 


R 


/ 


R 


J 


R 


I 


0.3 
0.5 
0.95 


Ord 
Ext 
Ord 
Ext 
Ord 
Ext 


-0. 0065 
-.0080 
-.0121 
-.0161 
-.2479 

-1. 1723 


-0. 9030 
1. 1074 

-0. 8676 
1. 1524 

-0. 3873 
1.8316 


-0. 0066 
-.0081 
-. 0123 
-. 0163 
-. 2505 

-1. 1928 


-0. 9027 
1. 1078 

-0. 8672 
1. 1530 

-0. 3838 
1. 8271 


0. 8859 
.6582 
.7983 
.2548 
.4634 
.9995 


-0. 0046 
-. 0439 
-.0091 
-. 2170 
-. 1160 
-. 7499 


0. 8858 
.6568 
.7982 
.2525 
.4625 

1. 0002 


-0. 0046 
-.0447 
-. 0091 
—.2224 
-. 1175 
-.7380 
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Next the values Y—l/2, Z m =0.l2 were used. These give ^ e ff~4.5X10 6 sec -1 which would 
apply at a height of about 70 km. The results are plotted in figure 4, and there is now a small 
but significant difference between the two formulas, which is greater for the extraordinary 
wave. This confirms Belrose and Burke's [1964] conclusion that for some purposes, at these 
low heights, the Appleton-Hartree formula may not be accurate enough. 

In figures 5 and 6 the values used are Y=6, Z m =l. Thus the curves refer to a frequency 
of about 200 kc/s and a height of about 70 km. In figures 7 and 8 the values are F=S0, Z m = 12 
so that they refer to a frequency of about 16 kc/s and again to a height of about 70 km. Fig- 
ures 5 and 7 give the refractive indices and show that in both cases the differences between 
the two formulas are significant though small. Figures 6 and 8 give the polarization p for the 
ordinary wave (sometimes called the "whistler mode"). The polarization for the other wave, 
according to (79), is simply the reciprocal of this. It is clear that the two formulas now show 
differences which are greater at the smaller values of X, that is below the level where these 
waves are reflected. 
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Figure 5. Similar curves to figure 4 but for the values 
Y = 6, e = 80°, Z m =l, Z ttf{ =2.S. 

These would apply for a frequency of about 200 kc/s and a height of 
about 70 km. 



Figure 4. A comparison of the values of the real parts 
At and imaginary parts — x of the refractive indices 
n according to the formula of Sen and Wyller (solid 
curves) and the Appleton-Hartree formula (broken 
curves) . 

In this example Y=H, 9=30°. For the solid curves Z m =0.12 and 
for the broken curves Zett =0.3=2.5 Z m . These values would apply 
for a frequency of about 2.5 Mc/s and a height of about 70 km. 
The difference between the n values for the ordinary wave is too 
small to show clearly on the scale of this diagram. 



O 

3- 0.5 




MINUS REAL PART 



12 3 4 5 6 7 

X - 



Figure 6. A comparison of the values of the real parts 
and imaginary parts of the polarization p for the 
ordinary wave according to the formula of Sen and 
Wyller (solid curves) and the formula of standard 
magneto-ionic theory (broken curves). 

The values used are as in figure 5. 
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Figure 7. Similar curves to figures 4 and 5 but for 
the values Y=80, = 30°, Z m =12, Z^f = 30, which 
would apply for a frequency of about 16 kc/s and a 
height of about 70 km. 



„ 0.5 




Figure 8. Similar curves to figure 6 but with the values 
used in figure 7. 



For these very low frequencies we have F>>1. Inspection of (68) shows that the argu- 
ments of the G functions in ei and e 2 are then approximately ±wY= ^z2irf H /p m . This is inde- 
pendent of frequency and even in the lowest part of the ionosphere (70 km) it is not less than 
3 or 4. If we therefore assume that |(1± F)w|>>1, the integrand in (67) may be expanded 
by the binomial theorem, and integrated term by term (compare (54) and (58)). If terms of 
the series beyond the second are neglected, the result leads to the values of standard magneto- 
ionic theory for e l and e 2 , given by (69). Thus only in e 3 should we expect a significant difference 
between the standard and modified theories. Figures 6 and 8 show that this affects the wave 
polarizations at low levels. 

At these low frequencies a change of polarization with height gives coupling between the 
two characteristic waves [Budden, 1961, ch. 19], and since the wavelength is large there can be 
a significant change of p within one wavelength. The reflecting properties of the ionosphere 
must usually be computed by a "full wave" solution of the differential equations, so that it is 
difficult to assess immediately the effect of a change of p on the reflection coefficient, but figures 
6 and 8 suggest that there might be a difference in the magnitude of the coupling process, ac- 
cording to the two formulas. The work of Deeks (private communication) has shown, how- 
ever, that the resulting effect on the ionospheric reflection coefficients is negligible for steeply 
incident waves in temperate latitudes. 

For waves whose wave normal is perpendicular to the earth's magnetic field, the refractive 
index for the ordinary wave is e\ /2 and is independent of ei and e 2 . In this case the curves of 
figure 3 could be used to find n, and the difference between the standard and modified formulas 
is more serious. This suggests that significant differences between reflection coefficients calcu- 
lated by the two methods might possibly be found for east-west propagation of VLF waves at 
grazing incidence, or at all angles of incidence near the magnetic equator. 

10. Effect of Stronger Wave Fields 

The foregoing results were derived with the assumption that the wave field E and the 
perturbation j x of the Boltzmann function / were so small that squares and products could be 
neglected. This may be called the first order approximation and an important consequence 
of it was that the isotropic part of/, that is/ , the spherical harmonic component of zero order 
in velocity space, was unaltered when the field was switched on, and was given by the Max- 
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wellian distribution (3). The average electron energy associated with the perturbation j x is 
zero because j x varies harmonically in time. Thus to the first order, the total energy of the 
electrons is unaltered. Between encounters the electrons are accelerated by the field E and 
acquire energy from it, but the generation of this energy is governed by a nonlinear term which 
was neglected when deriving (5). The energy acquired by the electrons in a very weak field E 
is too small to affect the current density J in (45), so that it can be ignored when calculating 
refractive indices, but it does cause a change of the second order in the isotropic part / 00 of the 
distribution function. 

To get the second order approximation to the Boltzmann function/, we must include 
squares and products of the small quantities j x and E, but we may continue to neglect higher 
powers. Since' the extra terms are small, it is accurate enough to use the first order approxima- 
tion in them. The results of doing this will be indicated here only in brief outline. 

Equations (4) and (17) are combined to give 

/(▼, *) = /oo(t>, *)+/ii(v, + /i 2 (v, t)+ (84) 

where the terms on the right are spherical harmonics of order 0, 1, 2 ... in velocity space. 
This is substituted in the Boltzmann equation (1). For the reasons given in section 3, we 
may neglect the grad r term which would allow for spatial dispersion and for inl Homogeneities 
in the medium. Then 

|ff£ E grad v/ll +^ (v A B) grad T /=|/ (85) 

where the first order approximation f n is used in the small second term. From this equation 
we now select the terms which are isotropic in velocity space (spherical harmonics of order 
zero), and equate them to zero. The third term makes no contribution, as was shown in 
section 5, so that 

^+±E.grad T / n =^. (86) 

at m & J ot 

This equation is not linear, so that the complex number convention for harmonically varying 
quantities cannot immediately be used, although it can be reintroduced as follows. 
When using the first order approximation f n , (20) must be replaced by 

gradv/n=J# [ { V+(v . |^) v/v | exp (ia*)] (87) 

and E must be replaced by 

${Eexp (iwt)}. (88) 

The second term of (86) is the product of the two harmonically varying quantities (87) and 
(88) both with angular frequency o>, so that it contains a part with angular frequency 2co, and 
a part which is independent of time and given by: 

i#(E*.V)+i#{(E*.v)(v.g)A>}. (89) 

Thus j/oo must contain a part with angular frequency 2a> and a more slowly varying part. 

When the radio wave is first switched on, the function / 00 changes as energy is imparted 
to the electrons. But these in turn begin to impart energy to the gas molecules, and eventually 
a new state of equilibrium is reached. The time constant for this process is of the order of 
1 msec [Huxley and Ratcliffe, 1949] and is important in the study of wave interaction. We 
now assume that the radio wave has been present for a long time, so that the slow variation 
of /oo has ceased, and only the variation of frequency 2co remains. Now let the time average 
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value of/oo be J 00 . Then-^=0, and the terms in (86) which are independent of time give 

de/oo _ /90) 



i^(E*.V)+i^^(E*-v)(y.g)/r}. 



dt 2 

Here V is given by the first order solution (40) or (42), and (90) is to be solved for the new 
distribution function y 00 . It should be particularly noted that E is a complex vector, and in 
general represents an elliptically polarized wave [Budden, 1964b]. A similar equation to (90) 
has been studied by Sen and Wyller [1960, appendix] but they have taken (88) above to be 
E cos ut where E is a real vector, so that their results apply only to linearly polarized waves. 
Next the left-hand side of (90) must be expressed in terms of f 00 and its derivatives with 
respect to v. To do this it is necessary to make a detailed study of electron-molecule en- 
counters, and it is not now permissible to assume that the molecules are infinitely massive. 
This is because the electrons are imparting to the molecules some of the energy which they 

have gained from the radio wave. For — ^^ Sen and Wyller used a formula derived by Chap- 
man and Cowling [1958, sec. 18.71] for electrons in a constant electric field, and this really needs 
modifying to allow for the harmonic variation of the field with time. Moreover this treatment 
took no account of molecular rotation caused by the electrons' impact. As far as the author 
knows, no proper treatment of this phenomenon has yet been given. 

The assumption in section 4 that the molecules are infinitely massive would mean that 
the energy of an electron is unaltered by an encounter. It would then be impossible for the 
electrons to impart any energy to the gas molecules, so that / 00 could never reach a constant 
value. But the more exact treatment indicated here shows that there is, on the average, a 
small amount of energy imparted to a molecule at each encounter. There is then a continuous 
flow of energy from the wave to the electrons and thence to the gas molecules. 

11. Conclusions 

This tutorial paper presents a derivation of the modified formulas of magneto-ionic theory, 
allowing for the dependence of the electron collision frequency on the velocity. No attempt 
has been made to set out the historical development of the subject nor to give references to 
all previously published work, and for further information the original paper of Sen and Wyller 
[1960] should be consulted. 

It is concluded that the standard Appleton-Hartree formula can be used without modifi- 
cation for nearly all radio propagation problems in the ionosphere provided that the value 
used for the collision frequency is the effective value v elt given by (60). This is %v m and is § of 
the average collision frequency. It is this v ett which is measured in those experiments which 
study the absorption of radio waves. For those measurements which use wave interaction 
[for example Fejer, 1955; Barrington and Thrane 1962] a more detailed theory will be needed 
on the lines indicated near the end of section 10, before any precise interpretation is possible. 
The effect of the modified formulas on the theory of the reflection of very low-frequency radio 
waves (10 kc/s to 100 kc/s) from the lowest ionosphere also needs to be studied for those cases 
where the wave normal is approximately perpendicular to the earth's magnetic field. 



The author is indebted to J. S. Belrose and B. H. Briggs for discussions which led to the 
writing of this paper, and to D. G. Deeks of the Radio Research Station, Slough, England for 
details of his VLF calculations. 
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